
STOP 



Early Journal Content on JSTOR, Free to Anyone in the World 

This article is one of nearly 500,000 scholarly works digitized and made freely available to everyone in 
the world by JSTOR. 

Known as the Early Journal Content, this set of works include research articles, news, letters, and other 
writings published in more than 200 of the oldest leading academic journals. The works date from the 
mid-seventeenth to the early twentieth centuries. 

We encourage people to read and share the Early Journal Content openly and to tell others that this 
resource exists. People may post this content online or redistribute in any way for non-commercial 
purposes. 

Read more about Early Journal Content at http://about.jstor.org/participate-jstor/individuals/early- 
journal-content . 



JSTOR is a digital library of academic journals, books, and primary source objects. JSTOR helps people 
discover, use, and build upon a wide range of content through a powerful research and teaching 
platform, and preserves this content for future generations. JSTOR is part of ITHAKA, a not-for-profit 
organization that also includes Ithaka S+R and Portico. For more information about JSTOR, please 
contact support@jstor.org. 



—67— 
ON AN UNSYMMETRICAL PROBABILITY CURVE. 

[Second paper.] 



BY E. L. DE FOBEST. 

In my article on this subject ending at p. 7 of Vol. X of the Analyst, 
some points of interest were omitted or but slightly touched upon. It is 
proposed to consider them here, and the formulas and tables obtained will 
be numbered consecutively with the previous ones. 

The accompanying figure has been drawn to represent the curve as found 
for the example at page 1, its ordinates being proportional to those in col- 
umn 9 of the table. For the equation of any curve of this nature, putting 
Ydx = y, we have bv (52), 

where a and b are constants determined by (39) from the q. m. error s and 
the c. m. inequality f, and K is a given function of a?b. 




The curve is in general limited on one side, where it touches the X axis 
at A, and extends to infinity on the other side, the axis being an asymptote 
to it in the direction OX. We may say, therefore, that it has a long side and 
a short side. Since in (81) we have 

x ax 

~al>~ a^b' 
and a s 6 is always positive, it appears that Y is a function of ax, and that its 
value is the same when a and x are both negative, as when they are both 
positive. Hence, changing the sign of a reverses the curve. 
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Our figure is drawn for positive a, but if it were negative the curve would 
be of the same form, only turned over so that the long side would come on 
the left hand and the short side on the right. From the nature of the curve 
as shown in (41) and elsewhere, it is evident that if the whole area included 
between it and the X axis could be cut out of cardboard or other material 
of uniform thickness, it would balance on a knife edge set to coincide with 
the Y axis. The greater leverage of the long side compensates for the great- 
er area of the short side. The vertex of the curve does not lie in the Y axis. 
The maximum ordinate CV is separated from OF by the interval OC, and 
this is important, for it shows that when the law of error of a set of obser- 
vations is such that £ is not zero, the arithmetical mean is not the most 
probable value of the observed quantity. To get the most probable value 
of a single observation, we must apply to the arith. mean a correction A 
equal to OC, and by (22) and (39) its amount will be 

A = - I = =£. (82) 

a 2e 2 

For the observations in Table II. we have 

A = — I = — = — 0°.91, 

a 1.095 

and this added to the arith. mean 5°.20 gives 4°. 29 as the value of the am- 
plitude most likely to occur. The probabilities that an error or deviation 
which occurs will fall on the -f- or — side of the most probable value, are 
not the same. The part VCX of the area of the curve lying on the long 
side is greater than the part AGV lying on the short side, as will appear 
hereafter. An ordinate dividing the total area into halves, would be be- 
tween OF and CV. The probable error of the most probable value of an 
observation is greater toward the long side of the curve than toward the 
short side. The ordinate which divides the area A CV into halves is PE, 
and that which divides VCX into halves is RF. The probable error on 
the + side therefore is CF, and that on the — side is CE. 

As found at page 3 and 4, putting n = a?b and v = ax, the areas on the 
short and long side of any ordinate Y whose abscissa x is reckoned from A 
as an origin will be respectively 

S] 7dm = 7^( 1+ ^ + (^+-2) +&C -> 1 (83) 
f> o-ie-" L . n— 1 , (n— l)(n— 2) . . \ f V ' 

and their sum is unity. The probability P 2 that an error which occurs 
will fall on the long side from VC is represented by the area VCX, and its 
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value is found from the second integral by taking as the lower limit x — ab 
— (1-7-a), whence we have 

v = ax — a 2 6 — 1 = n — 1, 

and the series terminates if n is an integer. The values of P 2 have been 
computed thus, to 7 places of decimals, for n = 4, 5, 6, 7, 8, 10, 14, 20, 30, 
50, 100, 200. The corresponding values of P x — 1 — P 2 , together with 
intermediate values found by interpolation, are given to 4 places in the sub- 
joined Table III. By means of the column of differences d lt P x and conse- 
quently P 2 can be found for any given value of the argument, with suffic't 
accuracy for all practical purposes. Now to get the position of the points 

Table III. 
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i?and F which mark the limits of probable error, we must find a; from each 
of the conditions 

C X Ydx = iP 1 , PYd» = lP„ 

^ o " " X 

and by (83) this amounts to solving the numerical equations 
u , = (l + _L_ +. .?! + &c.)~ ^•p i e>-*=0, 



= (l + ^=l + ^- 1 ] 2 ( w -=g) + &c ,\_(»- 1 )'p 2 e^-"^=0. J 






(86) 
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P,= 



value is found from the second integral by taking as the lower limit x = ab 
— (1-r-a), whence we have 

v = ax — a 2 6 — l = n — 1, 

- (^^^(2+^+M^+fe), (84) 
i (n) \ n — 1 (n — l) 2 / 

and the series terminates if n is an integer. The values of P 2 have been 
computed thus, to 7 places of decimals, for n = 4, 5, 6, 7, 8, 10, 14, 20, 30, 
50, 100, 200. The corresponding values of P 1 = 1 — P 2 , together with 
intermediate values found by interpolation, are given to 4 places in the sub- 
joined Table III. By means of the column of differences J lt P x and conse- 
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accuracy for all practical purposes. Now to get the position of the points 
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of the conditions 
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p 2 <?»— +1 = 



(85) 
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To illustrate the application of the table, take the amplitudes in Table 
II., for which we have obtained 4°. 29 as the most probable value. At p. 
3 we find a?b = 5.57. With this argument, Table III. gives 
= 1.085, ar 2 = 1.879; 

— ar, 1.085 



-ar-. 



- 1.095 
1.879 



— .99, 



= 1.72. 



a 1.095 

These are CE and GF in the figure. It is an even chance whether am- 
plitudes below 4°. 29 will be over or under 3°. 30, and an even chance 
whether those above 4°. 29 will be over or under 6°. 01. 

To show what appears to be the proper method of combining observa- 
tions when the c. m. inequality is of accidental origin, we will now find the 
most probable mean value and probable error of an angle observed with a 
theodolite by the method of repetition. The example is taken from Mey- 
er's Wahrscheinlichheitsrechnung, p. 274, where it is treated, of course, in 
the usual way, by the symmetrical law of probability. Fourteen observa- 
tions are taken, each giving the angle 17° 56', plus a number of seconds 
which stand in our Table IV. in the column headed u, The number of rep- 
etitions to each observation gives its weight, which stands in column p. 
The next column shows the products pu. The general arith. mean is 

1830 






_ ™v, = 39» 78. 
46 



Subtracting this from each u, we get the series of residuals v. 

Table IV. 
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to be reduced to unit weight, before combining them for the purpose of 
finding the law of error of w . According to (60), the q. m. error of unit 
weight is \/p times that of weight p, and the c. m. inequality of unit wei't 
is -fyp 2 times that of weight p. In other words, in a system of errors whose 
weig t is unity, the mean of the squares of the errors is p times, and the 
mean of the cubes of the errors taken with their signs is p 2 times what they 
would be in a system of errors having the weight p. Hence, as is well 
known, in order to find the approximate q. m. error e of unit weight from 
the squared residuals, each v 2 of weight p must be multiplied by p. Like- 
wise, to find the c. m. inequality £ of unit weight from the cubed residuals, 
each « 3 of weight p must be multiplied by p 2 . The last two columns in the 
table show the products pv 2 and pV. The number of observations being 
m = 14, we have for the square of the true q. m. error of unit weight 

; 2^l> 2 ] _ 1167.05 _ 8977 
1 m — 1 13 ' 

and as shown at p. 6, by virtue of the relation 

1 mm 2 ' 

the cube of the true c. m. inequality of unit weight is 

= m[ P V] = 14(-20396.6) = _ 14644> (89) 

1 m 2 — 1 195 

Then for the arith. mean, whose weight is [p] = 46, 

s: = ^= 89 f = 1.9515, ) 

C A ] _1 4 464 4 ^ ( 9 °) 

r* = ii_ = 14t>4 - 4 = —.6921. 

u [pj 46 2 J 

These are the (q. m. e.) 2 of u and the (c. m. i.) 3 of its possible errors, and 
from them we can construct the unsymmetrical cnrve which represents ap- 
proximately the probabilities of all the possible errors of u . 
By (39) we have 

2e n 2 /1.9515\ , RAr . \ 

« - ^ = 2 (~692l) = - 5 - 64 °' I (91) 

b = s 2 = 1.9515. j 

Since a is negative, the curve has the reversed position. By (82) the 
most probable error of u here is 

J > - ~« - if- •■•"•- jk = "- 18 - (92 > 

This is more probable than the error zero, or any other. Hence, for the 
best value of the mean we must take not u alone, but 

u -f- J = 39".78 + 0".18 = 39".96. 
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To find the probable errors of this result, we enter Table III. with the 

argument 

a 2 6 = 1.9515(5.640)2 _ q 2 .1, 

and there obtain 

—ar 1 = 4.892, ar 2 = 5.691. 

The probable errors are therefore 

— cwn 4.892 



-a 5.640 



0".87, 



_ _ ar 2 __ 5.691 __-,// 0] 
* 2 - "T ~=5^40 _ - U1 ' 

and the most probable value of the mean is written 

17°56'39".96+5qJ. (93) 

We have this final result instead of 17° 56' 39 // .78±0.94 as obtained by 
the ordinary method. Table III. also gives the probabilities that the error 
of m will fall on the short or the long side of the curve from u -f J , 
meaning in this case error in excess or defect respectively. They are 
P 1 = .4663, P 2 = 1— P x = .5337. 

While (93) shows an even chance for the error to fall within the limits 

+0.87 and — 1.01, the chance of its being between and 0.87 is not the 

same as that of its being between and — 1.01. The probabilities of the 

two events are 

JP X = .2331, |P 2 = .2669. 

We may observe that when the law of probability is unsymmetrical, the 
ordinary rule that probable error varies inversely as the square root of the 
weight, will not hold good. The probable errors depend not only on e but 
on £, and these, as (60) shows, vary with the weight by different laws. An 
observation of weighty is supposed to represent the arith. mean ofp obser- 
vations chosen at random, out of a set of observations of unit weight suffi- 
ciently numerous to show the unsymmetrical law of error. 

If the c. m. inequality of a set of observations of equal weight, is purely 
accidental, the probable value of C" is roughly expressed in termrs of e 3 by 
the relation (80). When observations are of unequal weight, £ 3 and e 3 va- 
rying with the weight in different ratios, we may consider (80) as approxi- 
mately valid for the average weight of the observations, so that s and £ 
denote the q. m. error and c. m. inequality for the weight \_p^-v-m. Then 
according to (60) 

" = (&}>>■ f! -fe)^ 

Substituting these in (80) we get 
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(Ci ») = ±:«^V(16[p]). (94) 

This is the limit within which there is about an even chance for ^ 3 to 
fall if the inequality is fortuitous. In our present example it gives 

(£•) = ± ^i 5 - T /[15x46(89.77)3] = 1076.4. 

The actual value by (89) is £/ = — 1464.4. The excess is not very 
great, and taken in connection with what we know of the nature of theodo- 
lite observations, it does not justify the belief that there was any real want 
of symmetry in the law of error. 

Granting that the true law is symmetrical, the most probable value of 
the observed quantity is the ultimate arith. mean, that is, the mean of all 
the possible values, each taken a number of times proportional to the prob- 
ability of its occurrence. But this ultimate mean cannot be found without 
taking a number of observations which is very large indeed, or infinite. 
We have to infer its probable value from a moderate number, in the forego- 
ing case only m = 14. Owing to the small number of observations, the 
distribution of the residual errors is irregular and unsymmetrical, and the 
law of error of the mean, so far as it can be deduced from the observations, 
is an unsymmetrical law. Hence, the most probable value of the ultimate 
arith. mean is not the apparent mean u , but u Q plus the correction J due 
to the c. m. inequality £ . The resulting value of the mean and its prob- 
able errors as exemplified in (93) constitute, as it seems to me, the best 
determination we can make of the most probable value of the observed 
quantity. 



CORRESPONDENCE. 



Editor Analyst: 

I desire to make an additional remark in regard to the paradox in the 
Query by Professor Johnson. (See Analyst, No. 2, p. 44.) 

The principle which seems to be violated in the Ex. given is this: — If w 
= f(a, x),.. (1) and du-t-dx =f'(a, x), . . (2) and we assign to a such value 
as will reduce u to a constant or to zero, then the same value ought to reduce 
/'(a, a) to zero. If it does not, the paradox must arise from a misinterpreta- 
tion of the result, or of some of the singular forms that enter into the work. 

Now although sin oo and cos oo may each be regarded as zero, in the 
sense that any constant may be regarded as zero in comparison with another 
quantity which is infinitely greater, we cannot admit as an absolute value 
in violation of the principle sin 2 c\3-j-cos 2 oo = 1. 



